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Abstract. The main objective of this article is to study the dynamic phase 
transitions associated with the spatial- temporal oscillations of the BZ reac- 
tions, given by Field, Koros and Noyes, also referred as the Oregonator. Two 
criteria are derived to determine 1) existence of either multiple equilibria or 
spatiotemporal oscillations, and 2) the types of transitions. These criteria gives 
a complete characterization of the dynamic transitions of the BZ systems from 
the homogeneous states. The analysis is carried out using a dynamic transition 
theory developed recently by the authors, which has been successfully applied 
to a number of problems in science. 



1. Introduction 

In 1950's, in his experiments, B. P. Belousov discovered a spatial-temporal os- 
cillation phenomenon in the concentrations of intermediaries when citric acid was 
oxidized by acid bromate in the presence of a cerium ion catalyst pQ. It is also 
observed by [2] that organic acids and metal ions could be used as well in the 
reaction, leading to spatial-temporal oscillations. It has been considered nowadays 
that all of the chemical reactions giving rise to oscillations, and the actions of cat- 
alyst are termed as the Belousov-Zhabotinsky (BZ) reactions. BZ reactions are 
now one of a class of reactions that serve as a classical example of non-equilibrium 
thermodynamics, resulting in the establishment of a nonlinear chemical oscillator. 

The main objective of this article is to study the dynamic phase transitions 
associated with the spatial-temporal oscillations of the BZ reactions, given by Field, 
Koros and Noyes [2J. This BZ reaction consists of the following five irreversible 
steps: 
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where 7 is a stoichiometric factor, P and Q are products which do not join the 
reaction again, and 

X = HBr0 2 , Y = Bv~, Z = Ce 4+ , A = B = BrO^. 

It is the temporal oscillation of the curium ion ratio Ce(IV)/Ce(III) which, with a 
suitable indicator, is displayed by a color change, when the reagent is stirred. 

The technical method for the analysis is the dynamic transition theory developed 
recently by the authors [6] [3 [8] . The main philosophy of the dynamic transition 
theory is to search for the full set of transition states, giving a complete character- 
ization on stability and transition. The set of transition states is represented by a 
local attractor. Following this philosophy, the dynamic transition theory is devel- 
oped to identify the transition states and to classify them both dynamically and 
physically. One important ingredient of this theory is the introduction of a dynamic 
classification scheme of phase transitions. With this classification scheme, phase 
transitions are classified into three types: Type-I, Type-II and Type-Ill, which, in 
more mathematically intuitive terms, are called continuous, jump and mixed tran- 
sitions respectively. The dynamic transition theory is recently developed by the 
authors to identify the transition states and to classify them both dynamically and 
physically; see above references for details. The theory is motivated by phase transi- 
tion problems in nonlinear sciences. Namely, the mathematical theory is developed 
under close links to the physics, and in return the theory is applied to the physical 
problems, although more applications are yet to be explored. With this theory, 
many long standing phase transition problems are either solved or become more 
accessible, providing new insights to both theoretical and experimental studies for 
the underlying physical problems. 

With this method in our disposal, we derive in this article a characterization 
of dynamic transitions of the BZ reaction. In particular, the analysis in this ar- 
ticle shows that the BZ system always undergoes a dynamic transition either to 
multiple equilibria or to periodic solutions (oscillations), dictated by the sign of a 
nondimcnsional computable parameter 5q — 5±; see (|3.12[) and p. 161) . 

For the case of transitions to periodic solutions (spatiotemporal oscillations), the 
Type of transitions (Type-I and Type-II) are determined again by another com- 
putable nondimensional parameter. In the multiple equilibrium case, for general 
domains, the transition is always mixed (Type-Ill), while for rectangular domain, 
the transition is either continuous (Type-I) or jump (Type-II) based again on the 
sign of another nondimensional parameter. To demonstrate the applications, the 
derived characterization of dynamic transitions of the BZ system is then applied to 
a special example. 

It is worth mentioning that continuous (Type-I) transitions imply that the con- 
centrations will stay close to the basic homogeneous state, and the jump (Type-II) 
transition leads to more drastic changes in the concentrations. The mixed (Type- 
Ill) transition lead to two regions of initial concentrations corresponding to jump 
and continuous transitions respectively. In addition, both Type-II and Type-Ill 
transitions are accompanied with metastable states, and fluctuations between these 
metastable states; see [11] for the related concepts for binary systems. 

This article is organized as follows. Section 2 introduces the basic model, and 
Section 3 study the dynamic transitions of the BZ model, and Section 4 gives an 
application of the theory and main results. 



PHASE TRANSITIONS FOR BELOUSOV-ZHABOTINSKY REACTIONS 



3 



2. FlELD-KOROS-NOYES EQUATIONS 

The stirred case was considered by [3], who derived a system of ordinary differ- 
ential equations for this reaction. Here we consider the general cases. Let X, Y,Z 
be variable, and A, B be constants. The equations governing fll.l[) are given by 

du 

—— = (JiAui + kiau2 — k 2 U\U 2 + k 3 bui — 2k i u\ 1 
at 

du 

(2.1) — — = (t 2 Aii 2 - k x au 2 - k 2 uiu 2 + jk 5 u 3 , 
at 

— ^ = a 3 Au 3 + k 3 bm - k 5 u 3 , 
ot 

where ui,u 2 ,u 3 ,a,b represent the concentrations of X, Y, Z 7 A, and B, Oi (i = 
1,2,3) are the diffusivities of and kj (1 < j < 5) are the reaction coefficients as 
in (|1.1[) . The model (|2.ip is also called the Oregonator. 

The coefficients Gi and kj are functions of the temperature T. In fact, kj = 
k^e~ Ei l RT , Ei is the activation energy, and R is the Boltzmann constant. 

The dimensions of the relevant quantities are given by: 

kj : M"H _1 for 1 < j < 4, k 5 : t' 1 , 

a, b, Ui : M for 1 < i < 3, Oi : \H" X for 1 < i < 3 

where t is the time, M is the mole density, and I is the length. Then we introduce 
the following nondimensional variables: 



k\a . k 3 b 
k 2 

k\k 3 



ui = - — m 1; u 2 
k 2 



u 3 = —^abuL t = {kik 3 ab)- 1/2 t' , 

k 2 k 5 

k 3 b ^ 
k\a 

/?=^T, S = k 5 (k 1 k 3 ab) 1 /", 



x = Lx' 

2ki 
k 2 k 3 b 



^ = P{k2ab)^ for i = 1,2,3. 

Omitting the primes, we obtain the following nondimensional form of (|2.1[) : 

dui 2 
= ^iAux + a(ui +u 2 - uiu 2 - pu^, 

(2.2) ^ = \x 2 Au 2 + -(ju 3 - u 2 - uiu 2 ), 

at a 

du 3 ,, . 

— = /i 3 Am 3 + - u 3 ), 

where the unknown functions are iti>0(l<i<3), and the parameters are 
positive constants: 

Mi, M2, M3, a, P, 1,5 > 0. 

Let f2 C ffi™, representing the container, be a bounded domain. 

If there is an exchange of materials on the boundary dQ to maintain the level 
of concentrations of X, Y, Z, then the equations (|2.2p are supplemented with the 
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Dirichlet boundary condition 

(2.3) u = (tii, u>2, U3) = on <9f2. 

If there is no exchange of materials on the boundary, the equations are supple- 
mented with the Neumann boundary condition 

Fin 

(2.4) ^ = on dSl. 

On 

It is known in [12] that the following region 

(2.5) D = {(«i,u 2 , "3) G i 2 (^) 3 I < <m < a u 1 < i < 3} 
is invariant for (|2.2l) - (|2.3p . where a, (1 < i < 3) satisfy 

01 > max{l, /3 -1 }, a 2 > 703, 03 > a\. 

The invariant region for the Neumann boundary condition case is the same as 
that of the Dirichlet boundary condition. For convenience, here we give the follow- 
ing lemma, which was well known in |121 113] . This lemma shows that the model 
has a global attractor. 



Lemma 2.1. The region D given by \2. 5)) is also invariant for the problem \2. 2\) 
with \2.J$. In particular this problem possesses a global attractor in D. 



Proof. We only need to prove that D is invariant for (|2.2[) with (|2.4p . It suffices to 
show that 

du 

(2.6) — points inward at x € dD. 



In fact, by (|2.2[) . we see that 

au2 > as ui = 0, U2 > 0, M3 > 0, 

-U3 > as U2 = 0, 112 > 0, M3 > 0. 
a 

Sui > as U3 = 0, ui > 0, U2 > 0, 
a[ai(l — /3ai) + (1 — 01)1/2] < as Hi = a±, U2 > 0, 

-[(02 - 7 U 3) + fltatii] < as u 2 = a 2 , u\ > 0, < u 3 < a 3 , 

a. 

-6(0.3 — ui) < as U3 = 03, < tii < o-i, 



dm 
~dT 

du 2 _ 7 
~dT 
du3 

dt 
du\ 

~dT 

du 2 \ 
~dt 

du3 

dt 



where a\ > max{l,/3 1 },a2 > 703, and 03 > a\. These properties above show that 
(|2.6|) is satished. The lemma is proved. □ 



Let 



= {(xi,--- ,x m ) € K m | > 0, 1 < i < m}, 
A = (/ii,/i2,M3,Q!,7,<y) € Mi. 
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The we define the following function spaces: 
H = L 2 {tt)\ 

!H 2 (fl) 3 n H^(Q) 3 for boundary condition (|2~3|) . 

f du 1 

< u e _ff 2 (£l) 3 I — = on <9i7 > for boundary condition 

Define the operators L\ = A\ + B\ and G\ : Hi — > H by 
A\u = (^iA-ui,^ 2 AM2,A i 3A-u 3 ) , 

f 1 7 

B\u = aui + au2, «2 H — W3, <5iti — (5i*3 

\ a a 

2^ 1 



G(u, A) = —a(uiU2 + Pui), U1U2, 

V a 

Thus, the Field-Noyes equations (|2.2p. with either (|2.3p or (|2.4p. take the following 
operator form: 

(2.7) ^=L x u + G(u,X), 

(2.8) u(0) = ip 

where A = (/xi,/z 2 ,/i3,o;,7j^) G R + . 

We note that the solutions u — (1*1,7/2,^3) of (|2.2[) represent concentrations of 
chemical materials. Hence only the nonnegative functions Ui > (1 < i < 3) are 
chemically realistic. 

3. Phase Transitions for BZ Reactions 

3.1. The model and basic states. We consider case where there is no exchange 
of materials on the boundary. In this case, the model is supplemented with the 
Neumann boundary condition (|2.4[) . and the system (|2.2p admits two physically 
realistic constant steady state solutions: 

(3.1) U = (0,0,0), I/i = («?,«§,«§), 



where 



u° 2 = ^ = hi +1 -M, 

1 + a A 



(1 - 7 - /?) + v/(l - 7 -/3)2 +4/3(1 + 7) 



It is easy to check that the steady state solution Uq = is always unstable, because 
the linearized equations of (|2.2p have three real eigenvalues Ai, A2 and A3 satisfying 




A!A2A 3 = det -i ^ =5(7 + 1) > 0, 



and hence there always at least one eigenvalue with a positive real part. 

Therefore we only need to consider the transition of (|2.2I) with (|2.4p at the steady 
state solution U\ in (|3.ip . For this purpose, we take the translation 

(3.2) u i = u' i + u° i (1 < i < 3). 
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Omitting the primes, the problem (|2.2j) with (|2.4|) becomes 



(3.3) 



-—- = fixAux + a[(l - u 2 - 2f3vfl)ui + (1 - u x )u2 - ^1^2 - /3u x ], 

5x^2 a ^ I" / 0\ 1 

— — = p 2 Au 2 + - \ -u 2 ui - (1 + + 7«3 - ""1"2 , 

at a 

du 3 

— = p 3 Au 3 + d(ui - u 3 ), 
du 



Then it suffices to study the phase transition of 



at u = 0. 



3.2. Linear theory and principle of exchange of stabilities (PES). In view 
of (|3.ip . the linearized eigenvalue equations of (|3.3[) are given by 



HiAui — a 



(7 + 3^(7- 1)U1 + ((7-1)1(2 



= Aui, 



(3.4) 



/i2Au 2 



1 



(1 + 7 — (3a)ui + (a + l)u 2 - "fu 3 



Xu 2 , 



p 3 Au 3 + S(ui - u 3 ) = \u 3: 
du , 



On 



an = 0. 



Let pk and be the fcth eigenvalue and eigenvector of the Laplace operator A 
with the Neumann boundary condition: 



Ae k = -ptek, {Pk > 0), 



(3.5) 

Let M k be the matrix given by 
(3.6) M k = 



de k 



Mi/»fc-f(7 + 3/3(7-1) 
-i(l + 7 -^) 

S 



-a{o - 1) 

"M2Pfc - + 1) 






J. 

a 

-P3Pk 



Thus, all eigenvalues A = (3 k j of Q3.4p satisfy 

M k x k j = PkjXkj, 1 < 3 < 3, fc = 1, 2, • • • , 

where € K 3 is the eigenvector of M k corresponding to /3jy. Hence, the eigen- 
vector u k j of (|3.4p corresponding to /3jy is 

(3.7) u kj {x) = x kj e k (x), 

where e k (x) is as in (|3 . 5[) . In particular, p\ = and e\ is a constant, and 



(3.8) 



-§( 7 + 3/3ct-1) -a(o--l) 

Mi= | -^(1+7-^) + I 

8 -S 
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£(3/J<7 + 7 -l) + -(a + l) 
2 a 



f (3/3(7 + 7 -l) + -(a + l) 



The eigenvalues A = /3ij (1 < j < 3) of Mi satisfy 
A 3 + AA 2 + BA + C = 0, 
A = S + 

(3.9) 

B = 2/3ct 2 + 7 - (1 - P)cr + (5 
C = 5a{2f3a + P + "f-l). 

It is known that all solutions of (|3.9[) have negative real parts if and only if 

(3.10) A>0, C>0, AB-C>0. 

Direct calculation shows that these two parameters A and C in (|3.9p are positive: 

(3.11) A>0, C>0; 
see also [4]. Note that 

PufaPia = -C<0. 

It implies that all real eigenvalues of (|3.8[) do not change their sign, and at least 
one of these real eigenvalues is negative. 

In addition, we derive, from AB — C = 0, the critical number 



(3.12) 
where 



c + b-a 2 + v /(c + 6-a 2 ) 2 +4a 2 & 



2 a 



a=£(3j9(r + 7-l) + i(<7 + l) 
2 a 

6 =(1-/3)^-2^2 -7, 
c= (7(2/3(7 + /3 + j- 1). 

It is then clear that So > if and only if 

(3.13) b = (1 - /3)(7 - 2fia 2 - 7 > 0, 

and under the condition (|3.13[) 



(3.14) 



> if 6 > S , 
= if S = S , 
< if S < S . 



Now we need to check the other eigenvalues fikj with k > 2. By 
(A; > 2) satisfy 



A fe =/3 : 



(3.15) 



A 3 + A k Xi + B k X k + C k = 0, 
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where 



A k =A + (/ii + p 2 + na)pk, 

B k =B + (fii + p 2 )pk + ap 3 p k + 

+ (MlM2 + MlM3 + P>2Pz)Pk; 

C k =c + 

\ a 

'1 + CT 



a + 1 . . . 

Mi + tt^aH + 7 - 1)M2 Pfe 

a 2 



cr + 1 a , 

a 2 



/0fc 



a 



Ml + o ( 3 ^ CT + 7 - l)/"2 M3 



pi + (J,l(J,2P3Pk, 



where a, b are as in (|3.12l) . and A, B, C are as in 
By (|3.1[) it is easy to check that 

3/3(7 + 7- 1 > 0. 

With the condition (|3~T3| we introduce another critical number 

PzPkb - C 



(3.16) 



5i = max 



^Mi + f (3/?cr + 7 - l)/i 2 ) p k + p x p 2 p\ 



P3pk 



where b is as in (|3.12[) . 

Then the following lemma provides characterizes the PES for Q3.4 



Lemma 3.1. Let S and Si be the numbers given by IS.lfy) and hS.lb}) , and b given 
by h3.12\) . When b < 0, all eigenvalues of \3.$ always have negative real parts, and 
when b > the following assertions hold true: 

(1) Let So < S\ and ko > 2 the integer that Si in 13.19$ reaches its maximum 
at pk - Then (3k i is a real eigenvalue of Ij3.4\ l, and 

< if S > Si, 

= if S = Si, forp k = p ko 
> if 5 < Si, 

Re/3ij(Si) < 0, # (k, 1) with p k = p ko . 

(2) Let So > Si. Then /3n(<5) = Pi2{S) are a pair of complex eigenvalues of 
\3.4ty near 5 = So, and 

( <0 ifS>S Q , 
RefSu = Ref3 12 I -0 if 5 = 5 , 

{ >0 ifS<S 0l 
Re(3 k3 (S )<0, V(fc,i)^(l,l),(l,2). 

Proof. In (|3.9[) we see that A = a + 5, B = aS — b. By the direct calculation, we 
can see that 

(3.17) A k > 0, A k B k - C k > 0, Vfc > 2. 

As b < 0, by (|3.11l) - (|3.13l) . (|3.15[) and (I3.17[) . for all physically sound parameters 
Mij M2j M3j a ? /^i <5 > 0, the following relations hold true 

A, > 0, C % > 0, AiBi - C z > 0, Vi > 1. 

Hence, it follows that all eigenvalues /3 k j(k > 1, 1 < j < 3) of (|3.8|) have negative 
real parts. 
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(3.18) RePij(Si) < 0, Vl<j<3. 



In addition, it is clear that 



(3.19) 



>0 if<5><5i, 
C kl { =0 i{S = Si, 
<0 H5<5i, 
Ck > at S = Si for all k ^ k\. 



Thus Assertion (1) follows from (|3~T7)) - (|3~T9)) . 

As Sq > Si, by p,19p we know that Ck > at S = So for all k > 2. Since the real 
eigenvalues of (3\j (1 < j < 3) are negative, the condition p,14|) implies that there 
are a pair of complex eigenvalues fin — /3i2 cross the imaginary axis at S = Sq. 
Then Assertions (2) follows. The lemma is proved. □ 



3.3. Dynamic Phase Transitions. By Lemma I3TTT as Si < Sq, the problem (I3.3|) 
undergoes a dynamic transition to a periodic solution from S = So. To determine 
the types of transition, we introduce a parameter as follows: 



o 2 3 

(3.20) h = [ Tr ((2£> - 6 - a 1P 2 D 3 )F 3 - (2jD 3 D 4 + aD 5 D - 6)Fi) 

- TT (("7P 2 ^3 - 2^ 6 )F! + (2 7 £ 3 ,D 4 ) + aD 5 D 6 )F 3 ) 

+ -~\-{2 1P D 3 + apD 6 + ajpD 3 D 5 - 2p- 1 D 4 D 6 )F 2 
F>o 



2D 2 Ep 



(27D3D8 + uDoD 7 )(2p- l DoD s - a 1P D 3 D 7 ) 



2 

{2 1 D 3 + aD 6 )(2 1 D 3 D 8 + aD 6 D 7 ) 



2D 2 E 

r,2 



^J2e(2D 6 - a 1P 2 D 3 ){2p- 2 DoD s - ai D 3 D 7 )] , 
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where 

p = Vs = y|, E = o?{a~lY(S 2 +p 2 ) 2 , 

n2 J 2 5 2 oP 2 (v I) 2 + P 2 ((S 2 + P 2 ) 2 + l8l{a l)) 2 j5 , a + 2<J 



iV 2 (a- 1) 2 (S 2 +p 2 ) 4 ' u aa 2 a(S-l) 

(3 (aA-2)(/3 ( 7 + 2/? + 7-l) n 2 - aA 

L>\ = 1 — — — , L>2 — 



a 2a(a-l) 2 ' ' a 2 (cj-l) 2 ' 

Z?4 = A + a/3cr 2 + a — 3a/3cr — a/? — 07, -D3 = So (a — 1), 

D s = A + a - a(3a - 2a/3 - ay, D 6 = (5q + p 2 ) 2 + 7^ 2 (ct - 1), 

Z? 7 = 1 - /3a - f3 - 7, D 8 = 13a 2 + (3 + 1 - 3/3a - 7, 
£>! 2p 2 J D 1 pD 2 

r\ — —. 



A A(A 2 + V) A 2 + Ap- 
„ D 2 Ap 2 D 2 2pD 1 

r 2 



2 ' 



^3 = 



A A(A 2 +Ap 2 ) A 2 +4p 2 
2pD 1 pD 2 



A(A 2 +Ap 2 ) A 2 +Ap 2 

Here A, B, C are as in (|3.9[) , and a is as in p.l2[) . 

Then we have the following dynamic transition theorem. 

Theorem 3.1. Let Si < Sq, and b\ is the number given by LS.20\) . Then the problem 
t3.3\) undergoes a transition to periodic solutions at S — Sq, and the following 
assertions hold true: 

(1) When bi < 0, the transition is of Type- 1, and the system bifurcates to a 
periodic solution on S < Sq which is an attractor. 

(2) When b\ > 0, the transition is of type-II, and the system bifurcates on 
5 > Sq to a periodic solution, which is a repeller. 

Proof. By Lemma 13.11 at 5 = Sq there is a pair of imaginary eigenvalues /3n = 
/3i2 = — ip of (|3.4|) . Let z = £ + in and z* = £* + in* be the eigenvectors and 
conjugate eigenvectors of (13.41) corresponding to — ip, i.e. z and z* satisfy that 

r . (M 1 +ip)z = 0, 

[6 > (Aff - ip)z* = 0, 

where Mi is the matrix as in (|3.8|) . and M* the transpose of M. Because ±ip are 
solutions of (13. 9L and AB = C at 8 — Sq, we deduce that 



(3.22) p 2 = B = Cj A. 

For z = {z\, z 2 , Z3), from the first equation of (|3.21[) we get 

(o 9 o\ S z 1 = (Sq - ip)z 3 , 

1 ' a(a - l)z 2 = (-§ (7 + 3/3a - 1) + ip) z x . 

Thus, we derive from Q3.23[) the eigenvectors z = £ + in as follows: 

(3 24) £ = (6,6,6) = (1, - ^Jf-i) 1 > 
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(3-26) ^ >~ yh l 



In the same fashion, we derive from the second equation of (|3.21|) we derive the 
conjugate eigenvectors z* = + irf as 

»?* = = (-^rry.0.- a(^+ P 2) ) ■ 

It is easy to show that 

= (»?.'?•) = -^. 

(£,*f) = -(r?, = 

where iii = a(cr — 1) (<5§ + p 2 ) 2 , D3 and Z?6 are as in (|3.20j) . It is known that 
functions ^P* +1^2 given by 

also satisfy the second equation of (|3.21|) with 



Let u e H be a solution of (|3.3[) expressed as 

it = + yry + y), x,y(zW}. 

where $(x, y) is the center manifold function of (|3.3|) at 5 = 5q. By (|3.28[) . the 
reduced equations of (|3.3|) on the center manifold are given by 



dx 1 

(3 29) i-« + ps ,G « + »' +ft ^ 

where G(u) = G(u, u) is the bilinear operator defined by 

(3.30) G(u, v) — (—auiV2 — affuivt, W1W2, 0) 

a 

for u = (ui,u 2 ,u z ),v = (ui,«2)W3) e fli. 

We are now in position to solve the center manifold function <fr(x,y). To this 
end, we need to determine the third eigenvalue /J13 and eigenvector £ of Mi at 
<5 = 5q. We know that 

As • (*p)(-ip) = V/?i3 = detAfi = C. 

By ([3~22]l we obtain 

/3i3 = --4 = -(5o + a), 
and a is the number as in (|3.12[) . Then, from the equation 

(M x - /? 13 )C = 0, 

we derive the eigenvector 

(3.31) C = (C, & ,C,)=(l /- f( ? + t + 1) .-^ 

\ a(cr — 1) a 

In the same token, from 

(Ml - p 13 )(* = 0, 
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we obtain the conjugate eigenvector as follows: 

(3.32) C* = (Ci.Ca.Ca) = ( < u — 

\ a(a — 1) aa 

On the other hand, from p. 241) and (I3.30|) it follows that 



(3.33) 



Gil = G(£,0 = (-«(&+ 0).-£>O), 
G 12 = G{t,ri) = {-am,-^ ,0), 
G 22 = G(r } ,r ) ) = 0, 
G 2 x =G(t},C = 0. 



Direct calculation shows that 

(C, C*) = A>, C) = D x , (G 12 ,C) = D 2 , 

and Do, Di,D 2 are as in p.20[) . 

By the approximation formula (A. 11) in |10j . the center manifold function <£> 
satisfy 



(3.34) 
where 



$ = $1 + $2 + $3 + o(2), 
£$1 = -x 2 P 2 G n - xyP 2 G 12 , 

(£ 2 + Ap 2 )C<5> 2 = 2p 2 (x 2 - y 2 )P 2 G n + Ap 2 xyP 2 G 12 , 
(C 2 + 4p 2 )$ 3 = p(y 2 - x 2 )P 2 G 12 + 2pxyP 2 G n , 

P 2 : H — > E 2 is the canonical projection, E 2 =the orthogonal complement of 
i?i=span{£, 77), and £ is the linearized operator of (|3.3p . By (I3.7[) . it is clear that 

(3.35) 



P 2 Gn = (Gu,C*)C = -DiC, 
P 2 G 12 = (G 12 ,C*)C = £>2C- 



Hence 4>i, $2, $3 £ span{£}. It implies that 



(3.36) 



£$, = Mi$ 7 



We infer from (|3.34[) - (|3.36[) the center manifold function as follows: 



(3.37) 



D 



A 



2p 2 D 1 



P D 2 



A(A 2 + 4p 2 ) A 2 + 4p 
>i 



4p 2 D 2 2pD x 



Do 



A A(A 2 +4p 2 ) A 2 +4p< 



( 2 P 



2 D, 



pD 2 



\A(A 2 +Ap 2 ) A 2 +Ap 2 



xy 
o(2) 



Do 



(Fix 2 + F 2 xy + F 3 y 2 )( + o(2), 



where Fi,F 2 ,F 3 are as in (|3.21[) . 



PHASE TRANSITIONS FOR BELOUSOV-ZHABOTINSKY REACTIONS 13 

Inserting ([337]) into ([3^9]) . by (|3~27|) . we have 

rl-r 1 

(3.38) — = - py+ _[[(£, £*)(G U ,0 + (£, ^(Gn, ?f )]x 2 

+ [(^n(Gi 2 ,n + «^*)(Gi2,??*)]^ 

+ ^-(£,O(G(£,0 + G(C,DO(^ 3 + ^ 2 y + iW) 
+ (£, r ) (G(C, v), C ) (^ 2 y + F 2 xy 2 + F 3 y 3 ) 
+ ^r(tv*)(G(L + G(C,0. ?7*)(^ 3 + ^ 2 y + iW) 
+ »?*)(G(C, »?), ^X^y + lW + F 32/ 3 )] + o(3), 

(3.39) ^ = P x+ i [[(»?, r)(Gii,r) + (»/,»l')(Gi 1) ») , )]i a 

+ [(r7,n(Gi2,r) + (r? ! r ? *)(G 12! r ? *)]^ 

+ ^-fa, D(G(e, C) + G(C, 0. 0(*i* 3 + fa^tf + F 3 xy 2 ) 
+ ^r(v, D(G(C, t?), D(^i^ 2 y + iW + F 3 y 3 ) 
+ ^-(??, ^)(G(C, C) + G(C, 0, ^)(^i^ 3 + F 2 x 2 y + F 3 xy 2 ) 
+ -£r(ri, V*)(G(t, V), V*)(Fix 2 y + F 2 xy 2 + F 3 y 3 )] + o(3), 

where D 2 = (CD 2 + (£,?f) 2 - 

Based on (ETM|) . (|3~251l . (|3"3D ]) -(|3"33 l) . we find 

( Gu ,0 = ^, (G 12 ,D = ^, 

X32 -"2 

(Gn,?? ) = 2J?2 , ( Gl2 > r l) = -^i 

(G(C,r ? ),n = -^, (G(C,r?),r ? *) = ^, 

J32 -"2 

(G(e, C) + G(C, 0,0 = ^. (G(C, C) + G(C, 0, - 

-fi 2 JJ 2 

Ha =a 2 (a-l) 2 . 
In view of (|3.26l) . equations (|3.38[) and (|3. 391) become 

■jT = - PV + [a 2 ox 2 + axxxy + a 30 x 3 + a 21 x 2 y + a 12 xy 2 + a 03 y 3 ] + o(3), 
^ =px + [b 20 x 2 + hixy + b 3Q x 3 + b 21 x 2 y + b 12 xy 2 + b 03 y 3 ] + o(3), 
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where 



p 2 

ZH1H2JJ 
a n = --^ r2 (2 7 D 3 + aD e ), 

620 = 2H 1 H 2 D 2 ( - 2aDeDs ~ Wp 2 D 3 D 7 ), 
611 = ^^(2^-^3), 

O30= D g 1 S 2 ^ (W4 + aD6Z?6) ' 

«12 = - ng p n2 [(27-P 3 J ) 4 + aD 5 D 6 )F 3 + (2 7 £ 3 + aZ) 6 )pF 2 ] 



DqHiH 2 D 2 
D Q H l H 2 D 2 



& 03 = n ; r ; r 3 n7 (2 As - aTp 2 /). 



3 J 



.,2 



6 2i = n »■ g na K 2 /'" D * D * ~ aipD 3 D 5 )F 2 + (2D 6 - orffi i D 3 )F L ] 

Due to Theorem A. 3 in [5], the number 

2 1 

h = 3(a 30 + 603) + (012 + &21) + -(002^02 - a 20 &2o) + -(ana 20 - &11&20) 

P P 

is the same as that given by (|3.20[) . Thus the proof is complete. □ 

3.4. Transition to steady state solutions. Thanks to Lemma T3.ll for <5i > Sq, 

the transition of p. 31) occurs at 8 = S\, which is from real eigenvalues. Let 5\ 
achieve it maximum at p ko (&o > 2). Assume that PkoliS) is simple near 8\. 

First, we consider the case where 



(3.40) f eljx^O. 

Jo. 



In general, the condition (|3.40j) holds true for the case where fi ^ (0, L) x D, and 
D C IR" _1 (1 < n < 3) is a bounded open set. 
Let 

1 a 
bo = [a/3(a - l)(p 2 p k(t + -(a + 1)) - (ap 2 p ko + 2<j)(p 1 p ko + —(7 + 3/3ct - 1))]. 

a 2 

Theorem 3.2. Let the above number bo 7^ 0, and Si > 8$. Under the condition 
the transition of i3. 3\) at 8 = Si is mixed (Type-Ill), and the system bifur- 
cates from (u, 8) — (0, Si) to a unique steady state solution u s such that u s is a 
saddle on 8 > 5%, and an attractor on 8 < Si. The solution us can be expressed as 



(3.41) 



u 5 = Cp kol (8)Ze ko +o(\p kal ), 
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where £ = (£1,^2, £3) and constant C are given by 
£1 = a(<r - l)(ti3Pk +<5i), 

Ct 

£2 = -OlPfco + ^(7 + 3/3cr - l))(fl 3 p ko + S^, 

£ 3 = aSi(a - 1), 



C = 



[(M3Pfc + ^i) 2 (^iPfc n + MgPfcp + j(o- + 1) + f (7 + - 1)) - (o- - l)7<5i] / n e fco riir 



Proof. We apply Theorem A. 2 in [5] to prove this theorem. Let £ and £* € M 3 be 
the eigenvectors of M ko and M£ o corresponding to P ko i{5\) = 0> i- e - 

M feo £ = 0, M* a C = 0, 

where Mfe is the matrix (|3.6I) with k = k - It is easy to see that £ is as in (|3.41l) . 
and 

/g £l = ~{^2Pk + ~(<T + l))(M3Pfc + <*l), 

£* = a(n3p ko +5i)(a - 1), 
£3=7(^-1)- 



For the operator G in (|2.7|) . we can derive that 



' -(G(^e feo ),re fco ) = -^ 2 + (2), 



(£e fc0 ,ee fe0 ) v v ^ ^ G^ 

where C is as in f|3.41|) . Therefore, the theorem follows from Theorem A. 2 in [9]. □ 

In the following, we assume that 

/ e 3 ko dx = 0. 
Jn 

Define 



(3-43) 6 1 =^i^[-( Ailj0fco + |( 7 + 3/3 ( T-l))(aM2Pfco+2) / o^'fj.r 





n 



+ 2a 2 f3(cr - l){p 2 Pk + — — ) / 0i4 o dx 
+ a(a-l){afi 2 pk +2) / <p 2 e 2 ko dx\, 



where 



(£,£*) = a(l-ff)[(WP*o +^2/5fc + ^^ + ^(7 + 3/3a-l)) 

a 2 

x(M3Pfc +<5i) 2 -7^1(0-- 1)], 
= (0i , 03) satisfies 

(3-44) L0 = -G(£)e 2 , 
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and the operators L and G are defined by 



L<t> 



~( 7 + 3/3<7-l)<h + (a -1)0 2 
i(l+7-/3(r)//i + (c7 + l)02 



G(0 



, M3A03 + (5(01 - 3 ), 

a 2 (fi - l)(/i3Pfc + <u) 2 MiPfco + 7^(7 + 3/3cr - 1) - af3(a - 1) 

(a - l){fi 3 Pk + Si) 2 (mp ko + -(7 + 3^cr - 1)), 
0. 



Then we have the following theorem. 

Theorem 3.3. Let b\ 7^ be the number given by \3.43\) , and 8q < 8\. Then the 
transition of h3.3\) at S = Si is continuous (Type- 1) as b\ < 0, and is jump (Type-II) 
as bi > 0. Moreover, we have the following assertions: 

(1) When b\ > 0, this system bifurcates from (u,S) = (0, 5i ) to two steady state 
solutions u s + and u s _ on 6 > Si, which are saddles, and no bifurcation on 
8<5 t . 

(2) When b\ < 0, this system bifurcates to two steady state solutions u s + and 
u s _ on S < Si which are attractors, and no bifurcation on 6 > Si. 

(3) The bifurcated solutions u± can be expressed as 

1/2 



± 



In e l dx 



-MS) 



£e fco +o(|/3 fcoi | 1/2 ), 



+ Si) 2 bi 

where £ € K 3 is as in Jff.^_?| ), /3 ko i(5) as in Lemma \3.1\ 
Proof. We use Theorem A.l in [9] to verify this theorem. Let 

and ^(y) is the center manifold function. Let $ = y 2 <fr = o(2); then by the 
approximation formula of center manifolds (see (A. 10) in [TU]), <fi satisfies 



= -G(Ze k0 ) = -G(0e 2 k0 , 



which is the equation (|3.44[) . 
We see that 



(G{y£e ko + $U*efc ) = - / [a{y£ie ko + $i)(y&e fco + $ 2 )£i*e fco 

Jn 

+ ap(y£ie ko + ^>i) 2 Cie ko + -(y^e ko + $i)(y&e ko + <S>2)iW ]dx 

a 

=y 3 [-« 2 £i* + ZafttQ + / ^iet dx 

+ / h4 dx] + o(3) 

a Jn 



Hence, have 
1 



(G(^ +$),fe ko ) 



In e l dx 



(fi2pk +Si) 2 biy 3 +o(y 3 ). 
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Thus the theorem follows from Theorem A.l in |9]. □ 

3.5. Stirred case. Theorem 13.11 describes the spatial-temporal oscillation for the 
BZ reactions of in a non-stirred condition. If the reagent is stirred, the equa- 
tions (|3.3p are reduced to the following system of ordinary differential equations: 

= -a Q(7 + 3(3<t - l)ui + (a - l)u 2 - u x u 2 - [3u\ 

(3-45) ^ = ™ f ~(1 + 7 - P<r)ui + (cr+ l)u 2 - -fu 3 - uxu^j , 

du 3 \ 
— =6{u x -u 3 ). 

In this case, only the transition to periodic solutions can take place, which is stated 
in the following theorem. 

Theorem 3.4. Let So > be the number given by i3.12)) . and b be as in \3.20\) . 

Then the system \3.4-5ty undergoes a dynamic transition to periodic solutions at 
5 = 5q. Furthermore, the system bifurcates to an unstable periodic orbit on 5 > Sq 
as b > 0, and to a stable periodic orbit on S < do as b < 0. In addition, the 
bifurcated periodic solution can be expressed in the following form 

(3.46) u s = [-b- 1 Reputf)} 1 / 2 ^ cos pt + rjshxpt) + o(|i?e/?n| 1/2 ), 

where £, r/ are as in A3. 32]) and \3.3$ , and /3n(<5) is the first complex eigenvalue as 
described in Lemma \3.1\ 

Remark 3.1. Since the constant eigenvector space E\ = spanjwn, u\ 2 , U13} = M 3 
is invariant for (13.31) where Uxj are given by (I3.7[l ). the bifurcated periodic solution 
u\ in Theorem 13.11 has the same equations of (|3 .45[) on the center manifold near 
5 = So, given by 

dr 1 

- = Ref3 11 (5)x - P (5)y + jrT^(G(x^ + yrf), £*), 

(3.47) f 1 

f t = p{S)x + Re(3 n (6)x + 7— r) (G(x£ + yrf), v *), 

where are given by ([3735]) and (f3736|) . p(S) = Im /3 n (S). The bifurcated 

periodic solution of (|3.45|) is written as 

u 5 = x(t)£ + y(t)r) + o(\x\,\y\). 

In the polar coordinate system 

a; = rcos6>, y = rshi9, 

the solution (x(t), y(t)) of (|3.47l) can be expressed by 

x(t) = a(6) cos pt + o(\a\), 

y(t) = a{6)~ pt + o(\a\). 

The amplitude a(5) satisfies 

Ref3 u (5) + ba 2 (S) + o(a 2 ) = 0, 
a(S) -> as 5 -> S , a(5) > 0. 

Thereby, we get 

a(6) = [-^Re^nWl^+otlRe/Jnl 172 ). 
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Thus, we get the expression (|3.46[) . 

Remark 3.2. By Lemma 12.11 for all physically-sound parameters, each of the 
two systems (|3.3|) and (|3.45j) has a global attractor in the invariant region D in 
(|2.5|) . Hence, when b > 0, the bifurcated periodic solution us is a repeller, and 
its stable manifold divides D into two disjoint open sets D\ and D2, i.e., D = 
D% + D2, D\ fl D2 = 0, such that the equilibrium Ui = (itj, u®, u®) in (|3.1[) attracts 
Di, and there is another attractor A2 C D2 which attracts £>2- 



4. An EXAMPLE 

We begin with chemical parameters. In the chemical reaction (jl.l[) , the constants 
given in [5] are as follows; see [J]: 

ki = lMM^S' 1 , k 2 = 1.6 x 10 9 M^S' 1 , k 3 = 8 x 10 S M- 1 S~ 1 , 

fc 4 = 4 x 10 7 A/ _1 5" 1 , a = i = 6 x 10" 2 M. 

Both 7 and k§ are order one parameters, and here we take: 

7 = 1, k 5 = l-S-\ 

and a = b as a control parameter. The nondimensional parameters a,/3,6 and [ii 
(1 < £ < 3) are given by 



k 3 b^ 



k\a 
2k\k^a 



= 7.727 x 10, 
3.375 x 1CT 6 , 



(4.1) ^ fc 2 /c 3 & 
5 = /c 5 (fcifc 3 a6) 1/2 = 1.035 x 10 2 aM _1 , 

w = » r ^ 7? = 4 - 664 x io-^m-^- 1 , 

(kik 2 ab)* I l a 

for i = 1,2,3. 

4.1. Stirred case. In view of (|4.1[) . the numbers in (|3 . 1 2[) are as follows: 

(4.2) cr = 7xl0 2 , a = 9.74, 6 = 690.79, c = 8.21. 
Hence 

(4.3) 5 = 71.67. 

Now we need to compute the parameter b\ in (|3.20[) . At the critical value 60 = 
71.67, the numbers in (13.91) are 



(4.4) A = a + 6 = 81.41, B = a8 - b = 7.27, C = S c = 588.41. 

Then we have 



p = 1 /B = 2.7, £^4xl0 21 , A) = 1.3 x 10~ 2 , 

Fi ^ 8.3 x 10~ 9 , F 2 = -2.6 x 10~ 7 , F a = -8.4 x lO" 10 , 

L> 3 = 5xl0 4 , L>4 = 4xl0 2 , A; = 80, 

Z3r = 3 x 10 7 D 7 = -6 x 10~ 3 , L»* = 4. 
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Thus, by flUTJ and flU^-flUS}, the number h in (f3T20]t is given by 

(4.6) 61 S x [("PAs - 2p- 1 D 4 D 6 )F 2 ~ aD 5 D e (3Fi - F a )] 



P 2 
L> 2 £ 



a 2 DlD 7 D 8 ot>p 2 D\D 1 3.2 x 10 7 D 6 Z? 8 



E£>V 2£L> 2 p 2 £;L> 2 



P 2 q 3.2D 6 D 8 7 
^ x 4 x 10 3 ^-r^ x 10 7 



EDqD 2 E 2 D 4 

6a 2 D 2 D 8 3 6aVD| 3 

^2^4 xlU - 2£;2jD4 xlU . 

Hence, in the stirred situation, by Theorem 13.41 and (I4.3[) - (|4.6[) , as S < Sq = 71.67, 
the system (|3. 451) bifurcates from (it, S) = (0, Sq) to a stable periodic solution, i.e. 
the reaction system undergoes a temporal oscillation in the concentrations of 
X = HBr0 2 ,F = Br", and Z = Ce 4+ . 



4.2. Non-stirred case. We consider the problem (|3.3[) . In this case, there is an- 
other critical parameter Si defined by (|3.16p . When Sq < Si, the system undergoes 
a transition to multiple equilibria at S = S±, and when Sq > S±, the system has a 
transition to periodic solutions at S = Sq. 

By (|4.ip and (|4.2p . the number o\ is given by 

T 690.8/i 3 588.4 

5l = max ol 777^ 1 7771 777^ — ^ 1 2 ~ ^Pk 

p k [9-1/ii + 0.7^2 + [iiii2Pk (9.1/ii + 0.7/i 2 )p/c + MiM2Pfc 

It is clear that if 

(4-7) M3<Mi,M2, 
then in view of (|4.3p we have 

690.8^3 690.8 r 

Oi < < < On. 

~ 9.1/ii + 0.7^ 2 9.8 
If (I4.7P is not satisfied, let f2n be a given domain, and f2 be the extension or 
contraction of 51n from xq in the scale L (0 < L < 00) defined by 

(4.8) 2:0) = {(x — x )L\ x, Xo £ O , Xq is fixed}. 

In this case, the eigenvalues pk (k > 2) of (|3.5p can be expressed by 

Pk = Tt" with C/t > and Ck —> 00 as k — >• 00. 

It follows that there exists an L > such that o\ < <5n for any L < Lq. Thus, by 
Theorem 13. II and (I4.6p . we have the following conclusion. 



Physical Conclusion 4.1. Under the condition of J^. 7| ) or SI = Q,{L, xq) with 
L < L Q , the system \3.S\) bifurcates from (u,S) — (0, So) to a stable periodic solution 
on S < Sq, which implies that the reaction system il.l]) with the Neumann boundary 
condition undergoes a temporal oscillation for S < Sq. However, when S > Sq, this 
system is in the trivial state U± in iS.l]) . 

Phase transition AT S = Si. By the above discussion, when > p,i(i = 1, 2) 
and L > L is large enough, the condition <5 < Si may hold true. In the following, 
we assume that L > Lq is sufficiently large, and 

Ml = A*2, M3 = zp,i (z > 1). 
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Then Si becomes 

r 691z 589 1 

S 1 - max — — — — xz , {x = p x p k ). 

x>o [10 + x (10 + x)x 

It follows that 3 < x < 4 for z > 1. We take x = pipk — 3-5, then 

Si = 47.7z- 12.5 

Thereby, we deduce that 

8i > 5 = 71.8 as — = z > 1.8. 

Mi 

Now, we investigate the transition of (|3.3p at 5±. For simplicity, we consider the 
case where 

n = (0,L) c K 1 . 
In this case, the eigenvalues and eigenvectors of (|3.5[) are 

Pk = (k - l) 2 7r 2 /L 2 ; e k = cos(fc - l)irx/L. 

It is clear that 

4 = (ei +e,-)/2 with j = 2*o-l, 

e| dx = by fc > 2. 

n 

We need to compute the number b x in [pT43|) . By (|4~Tj) . (j472l) - ((4~3|) . and 

M3Pfco = ^MiPfeo = zp2P/c = 3.5z (z > 1.8), 
the vectors G(£) and in (I3.44p are given by 

G(0 « (1.6 x 10 7 (/i3Pfe o + *i) 2 , 3 x 10 3 (/i 3 p fco + <5i) 2 , 0), 
= tp e x + Vjej (j = 2fc - 1), 

M^ = -~<?(0, 



2 



where 





1.25 x 10~ 2 
-(papj + Si) 

The direct calculation shows that detMi = — 8.755i, detMj ~ -1.5x 10 4 z+4.9x 10 3 , 
and 





! , 12.5*1 -5.4 x 10 4 S 1 * 
Mf 1 = r-rr 0.7*1 * 



detMi 



14/ 



! , 26.5(/i 3 pj +*i) -5.4 x 10 4 (Ai 3 Pj + *i) * 
M^ 1 = -1.25 x 10- 2 /i 3Pj 14.7( M3Pj + *0 
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Then we get 

<A) = (f 01,f02, V03) 
^02 = 1-1 X 10 2 (^ 3 Pfc o +<5l) 2 

13.4 x 10 3 



<p ai = 2.5xl0 6 (ii 3Pko +5 1 ) 2 



<Pji = 
Namely 



1.5z- 0.5 



<fj = {<Pjl,(Pj2,iPj3) 

2z + 27 



f]2 



1.52- 0.5 



01 = (poiei + = (^3Pfeo + <5i) 2 

02 = <A)2ei + ^ 2 ej = (fjiaPko + ^i) 2 



2.5 x 10% 



13.4 x 10 J 
1.521-0.5' 



Putting 0i and 02 into (|3.43[) , we derive that 

bi= (»s Pk0 +6ifL ( 5xlQS 5.7 xlO* 



, 2z + 27 

L1Xl0ei+ L5^^ 

2z + 27 



(£,£*) V 1.52-0.5 1.52-0.5, 

= -699a [16.7 x (S-l + 3.5z) 2 - 699<Si] . 

By Si > S = 71.8 and z > 1.8 we see that b\ < 0. Hence, by Theorem 13.31 the 
system fll.l[) undergoes a dynamic transition at 5 = Si to the following steady state 
solutions 

"* - n± (- 555^^ flwW )' / ' t " at T^ + '(^) ' 

which are attractors. 
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